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fii  B  S  T  R  A  C  T 

Acomputational  procedure  for  the  gravimetric  determination 
of  the  deflection  Of  the  vertical  is  presented  in  this  report. 
This  procedure  uses  the  square  development  exclusively. 
The  effect  is  computed  for  the  area  extending  from  the 
computation  point  to  the  outermost  limits  for  which  data 
is  available. 

'i'i 

The  mathematical  basis  for  the  computations  is  the  Stokes 
formula  for  the  undulation'Of  the  gWid,  and  the  companion 
formulas  of  Vening  Meinesz  for  the  deflection;  components. 
Deflection  formulas  of  Vening  Meinesz  are  given  in  the 
appendix. 

The  deflection  components  for  two  areas  were  determined 
using  an  electronic  computer  and  the  new  procedure.  The 
results  of  these  computations  were  compared  with  hand 
computations  which  were  completed  using  templates  with 
circles  and  radial  lines. 


FOREWORD 

>)  Improved  design  in  gravity  measuring  equipment  has^  resulted  In  the 
accumulatioh  of  gravity  data  in  an,  ever  ^nereasing  volume ,  Technology 
for  processing  such  data,  has  not  kept  pace  with  the  continually  improving 
means  of  data  gatherings.  In  the  past  gravity  data  has  been  processed 
manually,  by  a  technique  that  has  become  known  as  the  classic  solution, 
whereby  templates  consisting  of  circular  rings  and  radial  lines  are  used. 
The  accuracy  of  the  manual  solution  has  bqen  good,  but  the  method  is  so 
laborious  that  the  solution  of  a  great  many  points  would  require  an  excessive 
amount  of  time .  The  comi^tational  procedure  presented  in  this  report 
enables  the  rapid  processing  of  gravity  data  by  an  automated  technique.  The 
storage  requirements  of  an  electronic  computer  allow  and  encourage  an  in¬ 
crease  in  the  accumulation  of  gpravity  data .  The  computational  technique 
ensures  that  the  data  will  be  Used  for  the  determination  bf  the  deflection 
components . 
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PREFAGE 

The  need  for  the  computation  of  deflection  contponeiits  by  an  automated 
process  is  universally  recognized.  The  best  method  by  which  the  deflectioh 
components  can  be  obtained  is  perhaps  subject  to  debate.  The  method  pre¬ 
sented  in  this  report  attempts  to  standardize  storage  requirements  and  to 
minimize  the  need  for  total  number  Of  stored  gravity  values .  Accuracy  of 
the  square  method  of  computation  appears  to  equal  that  of  the  blassieal  ring 
procedure .  The  area  immediately  surrounding  the  iwiht  is  evaluated  pre  ¬ 
cisely,  rather  than  approximated.  The  zone  method  of  development  is.  Suited 
advantageously  to  solution  by  electronic  computer  . 

Appreciation  is  expressed  to  the  following  Individuals  who  have  contributed 
to  this  technical  report:  Mr  .  Abraham  Balak  for  his  valuable  technical  advice, 
Mr.  Nathan  Fishel  for  his  contribution  Of  comparative  coniputations,  Mr.  Albert 
McCahan  for  his  contribution  of  gravity  data,  and;  Mr  .  A.  E;  Craig  for  the  en¬ 
couragement  and  support  necessary  to  make  the  solution  possible. 
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A.  INTRODUCTION 

When  the  position  of  a  point  on  the  earth's  surface  is  detenhineci  by 
astronomic  observations,  the  results  are  referred  to  the  direction  of  gravity 
as  defined  by  the  leveling  device  on  the  observing  instrument.  Due  to  the 
irregular  distribution  of  masses  throii^out  the  earth,  the  level  surface  of 
reference  varies  from  place  to  place  in  an  irregular  manner,  The  gebdette 
computations  are  referred  to  the  regular  surface  Of  the  reference  ellipsoid.  ^ 

The  deflection  of  the  vertical  components  indicates  the  angular  tilt  between 
the  normal  of  the  ellipsoid  and  the  normal  of  the  level  surface . 

The  determination  of  the  deflection  components  has  always  be«;n  a  time-  ^ 
consuming  and  tedious  undertaking,.  This  report  describes  a  computational 
method  whereby  the  procedure  can  be  automated.  The  determination  of  the 
effect  includes  the  area  from  the  computation  point  to  the  outesr  limits  for  which 
data  is  available. 

0 

For  the  gravimetric  determination  of  the  deflection  components,  two 
general  methods  have  been  employed.  The  first  and  older  of  the  two  methods 
uses  templates  consisting  of  circular  rings  and  radial  lines  .  The  procedure 
consists  of  a  graphic  determination  of  the  anomaly  value  for  each  segpnent  df 
area.  The  anornaly  value  for  each  radial  segment  is  sunimed;  the  summation 
is  then  multiplied  by  a  factor  which  will  give  the  two  deflection  components. 

The  method  is  accurate  though  archaic,  since  each  step  has  to  be  performed 
manually.  The  method  cannot  be  automated  since  the  anomalies  are  for  area 
segments  that  have  a  unique  reference  to  a  given  computation  point. 

The  second  method  employs  a  technique  whereby  the  anomalies  are  per¬ 
manently  stored  for  areas  defined  as  geographic  quadrilaterals .  This  method 
is  referred  to  as  the  computation  with  squares .  The  greatest  advantage  Is  that 
the  anomaly  values  do  hot  have  to  be  re-estimated  for  each  deflection  computation. 
This  technique  is  well  suited  for  automation,  but  until  now  this  method  has  been 
used  only  for  the  determination  of  the  effect  of  the  outer  area .  The  circular 
method  has  been  used  to  complete  the  evaluation  of  the  inner  area.  The  disad¬ 
vantage  here  is  that  the  manual  circular  method  must  still  be  used  In  order  to 
finish  the  determination  of  the  deflection. 

The  method  presented  in  this  report  is  for  a  computational  procedure  accom  ¬ 
plished  by  the  use  of  squares  ^  The  main  disadvantage  of  previous  computation  by 
squares  has  been  overcome .  A  new  method  for  the  sumniatipn  of  the  anomalies 
for  the  inner  area  now  makes  it  possible  to  continue  the  deflection  computation 
rig^t  up  to  the  computation  point  •  The  advantage  is  that  the  whole  operation  can 
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be  eompletely  automated.  The  anomalies  are  permanently  stored  on  ma^etlc  tape, 
-  and  the  computations  are  performed  in  one  operation  on  the  eiectronic  computer. 

A  zone  system  of  squares  will  be  developed  for  the  determination  of  the 
deflection  effect  .  Every  square  In  a  given  zone  is  of  the  same  size;:  it  will  be 
shown  that  the  size  of  the  squares  in  any  given  zone  are  such  as  to  ensure  an 
accurate  deflection  determination. 

The  storage  of  anomalies  is  a  critical  consideration  for  this  type  of 
approach.  The  system  imposes  storage  requirements  which  must  be  satisfied 
in  order  to  Implement  the  new  method.  , 

A  new  method  must  be  compared  as  to  accuracy  and  economy  with  the 
method  that  is  presently  being  used.  The  hand  computation  method  which  ilses 
the  templates  will  satisfy  the  check  for  accuracy  and  also  provide  a  comparison 
with  the  new  method.  Two  problems  have  been  completed  by  both  methods  and 
the  results  are  presented  in  Section  G. 

>  o  ^  ■  •: 

structure  of  square  method 

A  zone  system  of  squares  is  developed,  whereby  72  squares  constitute  a 
zone,  and  each  square  in  a  given  zone  is  of  the  same  size.  See  Figures  1  and  2. 

A  zone, is  determined  by  dividing  a  given  square  area  into  81  smaller  squares; 
the  9  innermost  of  these  squares  are  deleted,  and  the  remaMng  72  squares 
constitute  the  zone .  The  next  zone,  approaching  the  center,  is  determined  by 
taking  the  9  squares  which  were  deleted  from  the  preceding  zone  and  subdividing 
this  area  into  81  squares  .  The  9  innermost  of  these  squares  are  deleted,  and  the 
remaining  72  squares  constitute  this  zone .  The  same  procedure  is  followed  with 
all  additional  squares.  ^ 

Zones  are  extended  inward  and  outward.  Any  square  in  a  given  zone  has 
the  relationship  to  any  square  in  an  adjacent  zone  such  that  the  sides  of  the  squares 
have  a  ratio  of  1:3  proceeding  Outward  and  3:1  proceeding  inward.  Starting  with 
Zone  1  where  the  72  squares  are  one  minute  squares,  and  proceeding  outward 
c  through  successive  zones,  a  total  of  8  zones  will  inelude  all  of  the  area  of  interest^ 
with  the  exception  of  the  inner  area.  Starting  with  Sub -zone  I,  which  is  just  inside 
Zone  1,  and  proceeding  . inward  through  successive  sub -zones,  an  infinite  number 
of  sub -zones  will  be  developed. 


''  o 


C  .  GENTER  OF  GRAVITY  OF  SQUARES 

The  squares  are  actuially  quadrilatera:ls  with  limits  defined  by  geographic 
coordinates  of  latitude  and  longitude.,  Figures  3a  and  3b  show  that  when  projected 
on  a  plane,  a  geographic  quadrilateral  can  be  represented  as  an  isosceles  trape¬ 
zoid,  Obviously ^  the  mean  latitude  and  mean  longitude  (Figure  3a)  will  not  locate 
the  center  of  gravity  of  a  quadrilateral , 


quadrilateral . 


In  Figure  3b  the  value  x  is  obtained  as  follows: 


(3)* 


and  the  value  of  y  is  determined  by: 

area  of  quadrilateral 

^  ,  d  +  X  '  (4) 


*  See  Appendix  for  expreSsldns  (l)„and  (2) 
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The  expressions  (3)  and  (4)  along  with  the  mean  longitude  define  the 
center  of  gravity  of  the  quadrilateral. 

If  the  mean  latitude  and  longitude  are  used  as  the  point  to  which  distance 
and  azimuth  are  determined  from  the  computation  point,  a  systematic  error  is 
introduced  into  the  computations.  To  eliminate  this  bias, the  center  of  gravity  of 
each  quadrilateral  is  determined  and  used,  The  following  tabulation  lists  the 
results  of  several  shifts  to  the  center  of  gravity  for  quadrilaterals . 


Example  I 

Example  2 

Emniple  3 

0, 

00°  00' 

o 

o 

o 

O 

CO 

60°  00’ 

01 

20°  00’ 

50°  00’ 

w.  w; 

M  0 

10°  00’ 

40°  00’ 

70°  00' 

0, 

09°  54’ 

39°’  21’ 

'67°  49' 

M  0  -  0, 

00°  06’ 

8 

Q 

02°  11’ 

D.  COMPUTATION  PROCEDURE 

For  Zones  1  through  8  the  computation  proceeds  zone  by  zone.  Each  zone 
is  completely  evaluated  before  going  on  to  the  next.  The  iimer  area  begins  with 
Sub -zone  1  and  continues  for  an  infinite  nuniber  of  sub -zones  approaching  the 
computation  point .  All  of  the  sub-zones  of  the  inner  area  are  evaluated  at  one 
time  by  the  summation  of  an  infinite  series . 

Figure  1  shows  that  the  72  squares  that  constitute  a  zone  have  12  different 
geometric  relationships  with  respect  to  the  computation  point.  These  12  relatlpn- 
ships  are  indicated  by  letters  a  through  1.  Table  1  is  given  as  a  proof  to  show  that 
regardless  of  sub -zone,  the  same  relative  square  in  each  sub -zone  exerts  a  con¬ 
stant  deflection  force  upon  the  computation  point. 


For  the  computation  of  the  inner  sub  -zones  a  constant  gradient  throughout 
the  area  is  assumed;  this  condition  is  more  nearly  true  for  each  successive  sub- 
zone  approaching  the  computation  point,  All  72  squares,  for  all  sub -zones, 
approach  the  central  point  in  precisely  the  same  mathernatical  manner.  Each 
square  approaches  the  point  at  a  ratio  of  two-thirds  with  respect  to  the  preceding 
square.  By  knowing  only  the  first  and  last  term  of  the  series,  the  summation  of: 
the  ahomalies  is  as  follows; 

5  Ag,  *=  Agj  +  Agj  +^,  ..Ag^  (5) 

where  successive  terms  of  Ag  are  determined  by  the  previous  md  last  term  aS  follows: 

Ag„  -  2/3  (  Ag,  -  Ag„.j)  +  Ag„^,  (6) 
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It  can  be  seen  from  (5)  that  if  Agp  is  anything  other  than  zeroi  then 
the  series  diverges,  and  the  sunimation  is  infinite.  Since  the  gradient  rather 
than  the  anomaly  Values  creates  the  deflection  effect,  the  whole  anomaly  pattefn 
of  the  iiiner  area  can  be  shifted  so  that  Ag^.  is  zero  •  This  can  always  he 
accomplished,  and  the  summation  can  always  be  rendered  finite  .  The  following 
expression  achieves  the  shift  and  summation  simultaneously  . 


OB 

SAg, 

lal  ' 


(Ag,  -  Ag„)  (1  -  0.333*) 
1  -  0.333 


(7) 


The  expression  (7)  reduces  to  (8): 


I^Ag,  =  1.5  (Ag,  - Agp  )  (8) 

where  Agj  is  the  anomaly  for  any  one  of  the  72  squares  of  Sub -zone  1,  and  Ag^^ 
is  the  anomaly  for  the  computation  point . 


With  the  summation  of  the  inner  area  achieved  by  (8),  it  is  only  necessary 
to  multiply  each  summation  by  the  constant  effect  . 

Table  1  was  developed  for  a  plane;  when  a  quadrilateral  is  used,  there  will 
be  72  constant  geometric  relationships  rather  than  the  12  shown  in  Table  1,  This 
does  not  impair  the  accuracy  nor  restrict  the  electronic  computer.  The  effect 
for  the  constant  geometric  relationship  was  tested  from  the  equator  to  latitude  8$* 
and  the  deviation  due  to  the  convergency  of  the  meridians  amounts  to  less  than 
1/100  of  one  percent. 


E.  EVALUATION,  OF  DEFLECTION  EFFECT  OF  A  SQUARE  ON  A  POINT 

The  area  shown  in  Figure  4  portrays  a  rapid  change  of  gradient .  If  the 
size  of  the  square: UVWX  Is  very  large  with  respect  to  its  distance  from  the  com¬ 
putation  point,  then  the  deflection  determination  fOr  the  square  will  be  in  error  by 
an  excessive  amount.  If  the  size  of  the  square  UVWX  is  very  small  with  respect 
to  its  distance  from  the  computation  point,  then  the  rapid  change  of  gradient  will 
be  insignificant,  and  the  deflection  determination  will  be  correct;  The  critical 
consideration  here  is  the  size  of  a  square  with  respect  to  its  distance  from  the 
point  .  If  the  squares  selected  are  too  large,  there  will  be  inaceuraCies .  If  the 
squares  selected  are  too  small,  there  will  be  too  many  of  them,  and  work  will  be 
increased.  It  is  the  purpose  of  this  section  to  show  that  the  size  of  the  squares 
used  are  such  that  accuracy  is  not  sacrificed. 


Figure  4.  Free-air  anomaly 
contours  showing  a  rajpid 
change  of  gradient. 

Assume  that  the  whole  area  of  square  UVWX  (Figure  4)  coincides  with 
square  52d  (Figure  1),  and  that  this  area  is  evaluated  directly  as  square  52d. 
Square  52d  is  then  divided  into  9  smaller  squares  and  each  of  these  Is  evaluated 
Independently  and  the  9  deflection  effects:  are  then  summed.  Square  UVWX  is 
then  successively  shifted  to  squares  48d,  2:ld,  and  25d,  In  each  of  these  loca¬ 
tions  (52,  48,  21,  and  25),  sc  uare  "d"  represents  the  same  geometric  relation¬ 
ship  with  respect  to  the  computation  point,  and  in  each  of  these  locations,  square 
UVWX  retains  the  same  orientation  with  U  always  located  at  the  upper  left- 
hand  corner  of  the  square .  The  same  computational  procedure  is  repeated  to 
square  UVWX  in  each  of  these  four  locations. 


Square 

52d 

48d 

2  Id 

25d 

Square  UVWX 

0V06G 

O'.'ObO 

0V060 

0VO6O 

9  Squares  X 

0VO62 

Q’.’Obl 

0V059 

OV06O 

A 

0'.'002 

OVOOl 

OVOOl 

OVOOO 

The  effect  of  square  UVWX  is  the  same  in  each  case,  since  a  mean  anomaly 
for  the  whole  square  was  used*  and  the  distance  and  first  quadrant  azimuth  ai^ie 
remaih  the  same .  The  effect  determined  by  the  9  smaller  squares  will  change, 
and  the  result  IS  a  more  precise  determination.  However  j  the  largest  difference 
for  A  is  only  QVO02,  and  it  should  be  noted  that  this  difference  is  found  in  an  area 
where  the  gradient  changes  rapidly  .  ; 


If  A  in  the  above  Gomparisons  had  been  large,  it  would  indicate  that  the 
square  "d"  was  too  large  with  respect  to  its:  distance  from  the  computation  point;. 
This  indicates  from  the  comparisdns  that  square  "d"  is  of  a  reasonable  size  arid 
that  good  results  can  be  obtained  even  when  a  rapid  change  of  gradient  dccufs 
within  the  square.  Every  square  in  every  zone,  in  the  method  presented  in. this 
report,  is  of  the  same  relative  size  with  respect  to  the  computation  point  .  The 
ratio  of  square  size  and  distance  to  the  central  point  is  a  constant^  Therefore, 
the  comparisons  given  would  be  valid  for  any  zone. 

F.  STORAGE  OF  ANOMALIES 

In  order  to  implenient  the  method  presented  in  this  report,  a  library  of 
stored  anomalies  must  be  established  which  is  compatible  with  the  computatiohal 
procedure .  It  is  envisioned  that  the  deflection  computations  will  be  for  either  the 
center  of  one-minute  squares  or  the  corners  of  one-minute  squares .  The  following 
storage  requirements  will  permit  the  deflection  computation  for  every  one-minute 
square  in  a  one-degree  block. 


Zone 

Area 

Storage  Square 

No.  of  Squares 

8 

180“ X  360“ 

5“  X  5“ 

2592 

6,7 

61*  x6l“ 

1“  X  1“ 

3721  , 

5 

23.6“  X  23.6“ 

10'  X  10' 

20164 

4 

10“  X  10“ 

5'  X  5’ 

14400 

1.2,3 

2.5“  X  2.5“ 

1'  X  1' 

22500 
t  63377 

The  total  number  of  squares  required  for  an  initial  effort  is  63377,  but  as 
the  library  of  stored  anomalies  is  increased,  the  requirements  are  more  easily 
satisfied.  For  example,  consider  a  oiie -degree  square  adjacent  to  another  pne  - 
degree  square  for  which  computations  have  already  been  performed.  The  additional 
number  of  quares  required  for  this  method  is  11250,  which  is  a  considerable 
economy  o/er  the  initial  requirements. 

G .  COMPARISON  OF  MACHINE  SOLUTION  WITH  CIRCULAR  TEMPLATE 
HAND  COMPUTATION 

Deflection  computations  for.two  specific  problems  were  computed  and  com  ^ 
pared.  The  first  problem  was  a  hypothetical  one  using  bathymetric  curves  as 
milligals .  The  second  problem  was  for  the  defleetton  effect  at  Columbus,  Ohio . 

The  gravity  data  used  was  a  free  air  anontialy  chart  that  had  been  prepared  at 
Ohio  State  University  , 
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Each  of  the  two  problems  was  computed  by  two  methods  -  the  new  square 
method  and  the  circular  template  method.  The  square  method  was  performed' 
by  ail  electronic  computer  (the  Bendix  G 15  -D),  and.  the  circular  method  was 
performed  by  hand  computation i  The  effect  of  the  outer  corners  iii  the  square 
method  was  deleted,  so  that  the  comparisons  could  be  niade  between  areas 
common  to  both  methods , 

In  the  hypothetical,  problem  (see  Figures  2  and  5),  the  outer  radius  for 
the  circular  method  is  2?025.  In  the  Columbus  problem,  the  Outer  radius  was- 
selected  aS  1?407  since  this  was  the  largest  radius  for  which  anomaly  values 
Were  available. 


Each  of  the  computations  by  the  circular  method  Was  performed  twice,  by 
two  individuals,  to  reduce  the  possibility  of  error  in  the  comparisons  ,  The 
results  of  the  comparisons  are  given  below. 

Hypothetical  Problem 

Method  Used  ^  ^ 


Circular 

GV345 

3’.’225 

Square 

0V321 

3V185 

A  0V024 

A  0V040 

Columbus  Problem 

Method  Used 

e 

n 

Circular 

-0’.’677 

-5V605 

Square 

-0'.'718 

-5'.'397 

A  -0V041 

A  -0,V208 

The  comparisons  between  the  automated  square  method  of  computation  and 
the  hand-computed  circular  method  indicate  that  accuracy  requirements  have  been 
met.  It  is  impossible  to  determine  which  method  is  more  accurate,  since  cUffer- 
ences  of  the  same  order  as  those  shown  will  Occur  when  two  individuals  using  the 
same  method  solve  a  given  problem, 
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H.  SUMMARY 

Ijumber  of  Computations .  The  machine  solution  procedure  is  hot  suited; 
advantageously  to  solving  single  deflection  computatioris;.  However,  machine 
solution  procedures  are  highly  desirable  when  large  numbers  of  deflections 
are  necessary. 

Accuracy  of  Data.  A  Gomparison  of  the  circular  and  square  coniputation 
methods  Indicates  that  the  results  of  either  method  are  of  the  same  order  of 
aeeuracy .  These  results  are  within  the  limits  of  absolute  accuracy  .  Absolute 
accuracy  for  the  deflectioh  determination  is  restricted  to  plus  or  minus  one  second 
of  arc .  The  principle  reason  for  this  restriction  is  that  large  areas  of  the  earth's 
surface  reniain  unsurveyed.  Figure  5  shows  that  the  circular -square  results 
are  comparable  because  the  individual  wedge  and  square  ^shaped  segments  are 
approxlniiately  the  same  size  .  Greater  accuracy  with  either  method  is  limited  by 
the  small  amount  of  complete  gravity  coverage. 

Use  for  Data.  Availability  of  large  amounts  of  data  would  permit  publi¬ 
cation  of  tables  and  graphic  contour  charts  for-  each  of  the  deflection  components. 

I.  CONCLUSIONS 

The  time  and  labor  involved,  in  computing  deflection  components  by  hand 
has  in  the  past  greatly  restricted  the  application  of  these  corrections  in  geodetic 
problems  .  In  the  case  of  triangulatiOh  alone,  the  deflection  components  can  be 
used  to  correct  the  observed  horizontal  angles ,  The  corrections  to  the  horizontal 
angles,  which  Can  be  as  large  as  several  seconds,  are  seldom  applied  because 
of  the  labor  involved  .  When  a  machine  solution  for  the  deflection  determination 
is  available,  many  refinements  can  be  made  in  many  different  kinds  of  work. 

The  machine  solution  necessitates  storage  requirements  which  are  reason^ 
able,  and  which  must  be  met  before  the  automated  method  can  be  achieved.  The 
constantly  increasing  volume  of  gravity  data  along  with  the  geodetic  requirements 
for  navigational  and  geodetic  satellites  makes  it  most  urgent  that  the  computation  „ 
be  automated. 


ST' 


13 


SELECTED  BIBLiOGRAPHY 

HEISKANEN,  W.  A;,  VENING  MEINESZ,  F.  A.  The  Earth  and  Its  Gravity 
Field,  McGraw-Hill  Book  Go.  ^  New  York,  1958  r 

HEl$KANEN,  W.  A.  ,  UOTILA,  U.  A.  Gravity  Survey  of  the  State  of  OMo, 
Diviaiion  of  Geological  Survey,  Department  of  Natural  Resduf 068^ 

State  of  Ohio. 

o  - 

HOSMER,  George  L.  Geodesy .  Second  Edltloh.  John  Wiley  aiid  Sons,  Inc., 
New  York,  1930:. 

KAULA,  W.  M. ,  "Accuracy  of  Gravimetrically  Computed  Deflection  of  the 
Vertical.",  Transactions .  Volume  38  Number  3,  pp  297-305, 

June  1957,  American  Geophysical  Union,  Washington,  D.C. 

QRLIN,  H. ,  'Deflection  of  the  Vertical  Above  Sea  Level,"  The  Journal . 

No.  7,  October  1957,  Coast  and  Geodetic  Suiwey,  U.  S.  Department 
of  Commerce,  Washington,  D.C. 

RICE,  D.  A.  ,  Deflection  of  the  Vertical  from  Gravity  Anomalies,  partial  copy 
of  report  to  C.  O.  AMS  1951.  ,  Coast  and  Geodetic  Survey,  U.  S. 
Department  of  Commerce,  Washington,  D.C, 

SCHREn'ER,  J.  B,,  Computation  of  Distances  and  Azimuths  with  Geocentric 
Elements,  Supplementary  Formulas  for  the  Coddington  Method., 
Technical  Paper  No .  36,  Mapping,  Charting  and  Reconnaissance 
Research  Laboratory,  Columbus,  Ohio,  1952. 

SOLLINS,  A.  D.  ,  'Tables  for  the  Computation  of  Deflections  of  the  Vertical 
from  Grayity  Anomalies Bulletin  Geodesique,  nouyelle  serie. 

No.  6,  pp  279-300,  1947,  Paris,  France. 


15 


APPENDIX  1 

F .  A .  Venlng  Melnesz's  Formulas  for  the  DefleGtion  Components 
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where  f  is  the  deflection  component  in  the  meridian 

is  the  deflection  component  in  the  prime  vertical 
G  is  the  mean  value  of  gravity  in  milligals  equal  to  979. • 

O  .  s> 

-  o 

■  :  Ag  is  the  anomaly  in  milligals 

^  is  the  angular  distance  from  the  computation  point 
a  is  the  geodetic  azimuth  from  the  south. 
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FORMULAS  FOR  THE  DETERMINATIONIDF  THE  UNDULATION  OF  THE  GEOID 

The  following  formulas  are  for  the  determinatioii  of  the  undulation  of  the 
geoid .  When  the  gradient  is  a  constant  in  the  area  about  the  eynSputation  point, 
the  undulation  effect  for  this  central  area  can  be  determined  with  formula  (15)  . 

The  Stoke  's  )Formula  for  geoid  undulation  is 
N  = 
where 

S  (^)  =  1+  csc-|  —  6sin-|^-5cos  i/r  —  3cosv»  lbggCsin^^  (10) 


For  a  given  compartment,  N  is 

N(<^.X)  =  ^  Ag(^.X)S(^)do-  (11) 

when  \j/  <  1°,  (12)  is  less  then  5%  in  error. 

'  (12) 

Letting 

N  is  ' 

N(,^.X)SCAg(,^.\)  ,  <14) 

The  summation  for  a  given  numbered  square  ,  for  successive  sub -zones 
approaching  the  computation  point,  may  be  achieved  with  (15) , 

O 

I  N  S  [agC+4,^(i)+  .  -  •  48„G(^”3  (15a) 

V  In  closed  form ,  ( I5ay  becomes 
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u  " 

"  |(^"®D‘8.+|Og,+  24g.)+-^3 

The  following  substitution  reduces  the  error  to  less  than  1%. 
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